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Abstract  -  A  quantitative  study  for  the  discrimination  of 
different  hepatic  lesions  is  presented  in  this  paper.  The  study  is 
based  on  the  fractal  analysis  of  CT  liver  images  in  order  to 
estimate  their  fractal  dimension  and  to  differentiate  normal 
liver  parenchyma  from  hepatocellular  carcinoma.  Four  fractal 
dimension  estimators  have  been  implemented  throughout  this 
work;  three  well-established  methods  and  a  novel 
implementation  of  a  method.  Analytically,  these  methods 
correspond  to  the  power  spectrum  method,  the  box  counting 
method,  the  morphological  fractal  estimator  and  the  novel 
modification  of  the  kth-nearest  neighbour  method.  The  Fuzzy 
C-Means  algorithm  is  finally  applied  revealing  that  the  k-th 
nearest  neighbour  method  outperforms  the  other  methods; 
thus  discriminating  up  to  93%  of  the  normal  parenchyma  and 
up  to  82%  of  the  hepatocellular  carcinoma,  correctly. 

Keywords  -  Hepatic  lesions,  fractal  dimension,  Brownian 
motion,  box-counting  method,  morphology,  k-nearest 
neighbour  method. 

I.  Introduction 

Improvements  of  different  medical  imaging  modalities, 
such  as  Computed  Tomography  (CT),  Magnetic  Resonance 
Imaging  (MRI)  and  Ultrasonography  have  dramatically 
increased  the  ability  to  detect  and  diagnose 
liver  abnormalities.  Nevertheless,  biopsy,  an  invasive 
technique,  remains  the  most  effective  for  the 
characterisation  of  different  liver  abnormalities.  To  limit  this 
process,  new  techniques  from  the  image  processing  field 
may  be  applied  to  liver  images  in  order  to  isolate  regions  of 
interest  and  detect  liver  abnormalities. 

Consequently,  ultrasound  B-scan  imaging  has  become 
one  of  the  most  popular  modality  to  image  human 
abdominal  organs  such  as  liver.  These  images  appear  as 
textural  ones  and  various  techniques  have  been  applied  in 
order  to  characterise  different  tissues  to  normal  or  abnormal 
[1],  Recently,  several  investigators  used  MRI  and  CT  in 
order  to  evaluate  local  hepatic  lesion  of  all  varieties  [2], 

Detection  of  hepatic  lesions,  using  the  aforementioned 
imaging  modalities,  has  been  realised  from  the  processing  of 
liver  images  using  different  imaging  processing  techniques, 
such  as  texture  analysis,  grey-scale,  shape  descriptors,  etc. 

[3] ,  In  particular,  texture  features  based  on  the  Fourier 
power  spectrum,  the  grey-level  different  statistics,  the  grey- 
level  run-length  statistics,  the  spatial  grey-level  dependent 
matrices  have  been  used  for  the  liver  tissue  characterisation 

[4] ,  Most  of  these  liver  tissue  differentiation  techniques 
were  based  on  the  analysis  of  liver  images  with  the 
emphasis  on  the  problem  of  classification  [4], 


Moreover,  fractal  geometry  has  received  much  attention 
as  a  useful  tool  for  image  analysis  [5],  The  intensity  surface 
of  an  image  can  be  considered  as  a  fractal  object  whose 
properties  are  quantified  numerically  by  the  use  of  the 
fractal  dimension.  For  an  image,  the  fractal  dimension  is  a 
non-integer  number  between  2  and  3  and  it  is  a  measure  of 
the  roughness  of  its  intensity  surface.  Experiments  have 
demonstrated  that  the  fractal  dimension  is  highly  correlated 
with  the  human  perception  of  image  texture;  the  rougher  the 
texture  appears  the  larger  is  the  fractal  dimension. 

The  present  paper  consists  of  a  comparative  study  of 
four  fractal  dimension  estimators  in  order  to  discriminate 
different  liver  lesions.  These  fractal  estimators  correspond  to 
three  well-known  techniques,  such  as  the  power  spectrum 
method,  the  box  counting,  and  the  morphological  fractal 
estimator  and  a  novel  modification  of  the  kth-nearest 
neighbour  method.  These  fractal  dimension  estimators  have 
been  implemented  and  applied  to  CT  liver  images  to 
differentiate  the  following  cases:  normal  parenchyma  and 
hepatocellular  carcinoma. 

II.  Methodology 

The  analysis  presented  in  this  study  was  performed  in 
three  main  steps:  (a)  Data  acquisition,  (b)  Image  pre¬ 
processing,  and  finally,  (c)  Estimation  of  the  fractal 
dimension. 

A.  Data  Acquisition 

All  CT  images  were  captured  using  a  Philips  LX  CT 
scanner  of  the  Department  of  Radiology,  at  the  Eugenidion 
Hospital  -  University  of  Athens,  Greece.  The  scanner 
consists  of  a  major  system  of  a  Local  Area  Network  (LAN) 
and  it  is  connected  to  an  Agfa  Impax  MC-300  Medical 
Gateway,  a  system  that  digitises  the  CT  images  throughout 
the  Video  Spot  Imaging  (VSI)  component.  As  a  result,  CT 
images  of  512  x  512  pixels  and  256  grey-level  distributions 
were  digitised  and  then  driven  to  a  Dicom  server,  according 
to  the  Agfa  Protocol  IP.  Finally,  the  images  were  sent  to  a 
PC  with  Windows  NT4.0  for  further  processing. 

In  the  present  study,  two  sets  of  liver  images  were 
finally  acquired  corresponding  to  99  normal  parenchymas 
and  50  hepatocellular  carcinomas. 
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B.  Image  Pre-processing 

The  selected  abdominal  CT  images  were  then  processed 
in  order  to  extract  a  particular  region  of  interest  (ROIs). 
These  ROIs  have  been  manually  drawing  by  a  team  of 
radiologists  so  as  to  contain  only  liver  parenchyma  (normal 
or  abnormal),  with  no  major  blood  vessels  information.  The 
ROIs  were  selected  along  the  centre  line  of  each  image  and 
were  consisted  by  rectangular  areas  from  25x25  up  to 
88x88  pixels  size,  depending  on  the  liver  image. 

There  were  cases  where  ROIs  are  difficult  to  be 
discriminated  from  the  surrounded  liver  tissues.  Therefore,  a 
contrast  enhancement  method,  such  as  the  histogram 
equalisation  technique  [6]  was  priory  performed,  in  order 
the  ROIs  to  be  successfully  identified  by  the  experts.  It  must 
be  pointed  out  that  the  calculation  of  the  fractal  dimension 
estimators  were  performed  on  the  original  data  and  the 
enhanced  process  had  been  used  only  for  visualisation 
purposes. 

C.  Fractal  Dimension  Estimators 

Several  methods  have  been  proposed  for  the  estimation 
of  the  fractal  dimension  of  liver  images.  In  the  present 
study,  the  fractal  dimension  of  the  various  ROIs  had  been 
computed  throughout  four  different  estimators,  regarded  as 
the  most  representative: 

•  The  Power  Spectrum  Method  (PSM), 

•  The  Box-Counting  Method  (BCM), 

•  The  Morphological  Fractal  Estimator  (MFE), 

•  The  kth-Nearest  Neighbour  estimator  (K-NN) 


The  Box-Counting  Method  (BCM) 

The  box-counting  method  estimates  the  fractal 
dimension  of  a  signal  (image)  as  an  upper  limit  of  the 
Hausdorff-Besicovich  dimension  [5].  The  box-counting 

dimension  of  a  set  S  <z  R"  is  defined  as  follows: 

FO.lim108^  (4) 

log(l  /  r) 

where  N(r)  denotes  the  number  of  n -dimensional  cubes,  size 
r,  needed  to  cover  set  S.  The  image  plane  (mm)  is  covered 
by  a  3 -dimensional  grid  of  cubes  for  various  grid  sizes  r. 
The  number  of  cubes,  N(r),  containing  at  least  one  pixel  of 
the  image  is  counted  and  the  fractal  dimension  is  obtained 
by  the  slope  of  the  best  fitting  line  at  the  points 
(-  logr.logjV(r)). 

The  Morphological  Fractal  Estimator  (MFE) 

The  morphological  fractal  estimator  belongs  to  the  area 
measurement  methods  [8].  The  covering  methods  can 
measure  the  Minkowski-Bouligand  fractal  dimension  of  a 
fractal  signal  by  creating  multiscale  covers  round  the 
signal’s  graph.  More  precisely  these  methods  are  based  on 
the  multiscale  morphological  erosion  and  dilations. 

The  dilation  and  erosion  of  a  function  f. 
Df=[ 0,xjx  [O,  X,]— >  91  by  a  structuring  element  ge: 

[)  _>  are  defined  as  follow: 

(/©  g£)<X>*2)  =  SUP{/<X  ~S1>X2  ~52)  +  g(51,52)/ 

(xx  sx,x2  s2)£  Dj-,(sx,s2)e  1 


The  Power  Spectrum  Method  (PSM) 

The  power  spectrum  method  belongs  to  the  fractional 
Brownian  motion  (fBm)  methods.  The  image  is  assumed  to 
be  fBm  [7]  with  parameter: 

H  =  3-  FD  (1) 

with  0  <  H  <  1 .  Then,  the  power  spectrum  density  of  image 
is  given  by: 


(2) 


where  k  is  a  positive  constant.  The  exponent  b  is  related 
with  the  fractal  dimension  as  follows: 


b  =  2  +  2H  =  2(4- ED)  (3) 


(f  ©  ge)(xl,x2)  =  h\f{f(xl  +  sx,x2  +s2)-g(s1,s2)l 

(x,  +sx,x2  +s2)e  Df,(sx,s2)&Dgt ) 

where  ©  denotes  dilation  and  ©  erosion  and  (x, ,  x2 )  e  D  f  . 

Let  f(mx,m2)  be  a  two-dimensional  signal  where 
(m,  =0,1,. ...Mj  —  l,m2  =  0,1,..., M,  -1) .  The  fractal 
dimension  is  then  estimated  as  follows: 

STEP  1:  Select  a  convex  set  structuring  element  Bd  to  be  a 
discrete  version  of  a  continuous  B  which  satisfies  the 
theorem  of  [8]  and  for  which  the  gE  structuring  elements  are 
constructed. 

STEP  2:  For  e  =  1,2,...,  £max  the  parameters 

Ue (ml,m2 )  =  (f  ®ge )(m{, m2 ) , 

Le  (mx ,  m2 )  =  (/  ©  ge  )(m{ ,  m2 ) 
are  recursively  calculated: 


where  2<  b  <  4  .  Pentland  [7]  estimated  the  exponent  b  for 
various  directions  of  the  Fourier  plane  as  the  slope  of  the 
least  squares  line  at  the  points  (-logf,logP(f!  ,f ,  )) .  These 

estimates  were  then  collapsed  into  one  average 
measurement,  from  which  the  fractal  dimension  was 
obtained. 


Ue= i  (m,,  m2 )  =  (/  ©  ge=1 )(»;, ,  m2 ) 

=  max  {f(m1  +i,m2  +  j)  +  ge=1(i,j)i 

Ue  (m, ,  m2 )  =  (t/£_!  ©  ge=l  )(ml ,  m2 ),  £  >  1 
and 
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4=i Oi > nh )  =  (/  ©  gE=i )Oh ,m2) 

=  min  {/(«!  +  i,/n2  +  y)  -  g£=1  (i,y)} 

— i</,y<i 

4  (/Wj ,  m,  )  =  (Z£_!  ©  g£=1  )(/m,  ,  m2 ),  £  >  1 

STEP  3:  Compute  the  parameter  vol(Ag  ("£))  according  to 
the  following  equation: 

vo/(t,(£))  =  ^'^Ue{ml,m2)-Le(jnvm2)  (5) 

m,  =0m2  =0 

STEP  4:  The  Minkowski  -  Bouligand  dimension  is  obtained 
by  the  relation: 

DMB=3-a  (6) 

where  a  is  the  slope  of  the  best  fitting  line  at  the  points 
(lne,ln  VO  l{Ag  (£))),  for  s=  1,2,...,  emax. 

The  kth-Nearest  Neighbour  estimator  (K-NN) 

In  a  previous  work  [9],  an  algorithm  for  calculating  the 
fractal  dimension  of  grey  level  images  using  the  kth-nearest 
neighbour  approach  was  presented.  In  this  work,  a  slightly 
modified  method  is  applied  for  the  estimation  of  the  fractal 
dimension. 

In  particular,  let  l(x,  y)  denotes  a  grey  level  image  with 
size  Nx  xNr  pixels,  then  the  fractal  dimension  is  estimated 
iteratively,  using  the  following  equation: 

<  r/  >~  (k  /  N )r/z)(7)  (7) 

for  k  =  kmtn £max  ( k  integer),  where  <  rk  > is  the  scaling 
of  the  average  distance  of  a  point  to  its  kth  nearest 
neighbour  as  a  function  of  k  ,  N  the  number  of  points,  and 
£>(y)  the  dimension  function.  The  fractal  dimension  is  then 
estimated  as  follows: 

STEP  1:  An  initial  value  of  y,  y0  =  2.5  ,  is  chosen 
arbitrarily. 

STEP  2:  Each  pixel  of  the  image  with  spatial  coordinates 
(jt,  y)  is  considered  as  a  point  of  R3  with  coordinates 
(x,  y,l{x ,  _>’)) .  For  each  such  point,  which  will  be  referred  as 
a  reference  point,  its  kmin  up  to  kmax  nearest  neighbours  are 
recorded  as  rk  (in  =  1,  2, . . . ,  N) ,  where  N  =  N x  xN 

denotes  the  total  number  of  the  pixels  in  the  image. 

STEP  3:  For  n  =  1, 2, . . .  the  following  recursive  relations: 

o(r>— ,7„=©(y„)  (8) 

Sn- 1 

are  applied  until  convergence  is  achieved,  s n_t  is  the  slope 
of  the  best  fitting  line  (using  linear  regression)  at  the  points 

(log(k/ N),  log  <  rl" >)  and  <  ryk" '  >= — ■  The 

distance  between  two  points  is  determined  using  the 
Euclidean  norm.  The  process  can  be  terminated  when  the 
fractional  range  between  ynl  and  D(yn), 


T>(7„)-y„- 1 

0-5[o(y„)+  y„_,] 


is  smaller  than  some  predefined 


tolerance  (e.g.,  10  5 ). 


Ill  Results 

The  data  used  in  the  present  study  are  composed  of 
rectangular  ROIs  of: 

•  Normal  parenchyma:  99  ROIs. 

•  Hepatocellular  carcinomas:  50  ROIs. 

In  Fig.  1,  two  typical  CT  liver  images  with  or  without 
hepatic  lesions  are  presented.  The  rectangular  areas 
correspond  to  the  best  fitting  rectangular,  within  the  ROIs, 
as  drawn  by  the  experts,  where  the  fractal  dimension  is 
estimated  using  the  aforementioned  methods. 


(a)  (b) 


Fig.  1.  Typical  transverse  CT  slices  of  (a)  normal  liver  parenchyma 
and  (b)  hepatocellular  carcinomas.  The  rectangular  areas  represent  the  ROIs 
under  study. 


In  Table  I,  the  mean  and  the  standard  deviation  values 
of  the  fractal  dimensions,  using  the  four  aforementioned 
fractal  dimension  estimators,  and  for  each  hepatic  case,  are 
then  presented. 


Table i 

THE  MEAN  AND  THE  STANDARD  DEVIATION  (IN  PARENTHESIS) 
OF  THE  FRACTAL  DIMENSIONS  FOR  ALL  ROIs  AND  FOR  EACH 
HEPATIC  LESION. 


PSM 

BCM 

MFE 

K-NN 

Normal 

Hepatocellular 

carcinomas 

2.25  (0.47)* 

2.85  (0.17)' 

2.18  (0.09) 

2.14  (0.15) 

2.86  (0.06) 

2.86  (0.06) 

2.26  (0.07) 

2.04(0.11) 

The  values  (mean,  standard  deviation)  of  the  fractal  dimension  for  the 
PSM  method  were  calculated  for  a  substantial  number  of  the  original  data 
(normal  parenchyma:  57  and  hepatic  carcinoma:  46)  since  the  method 
underestimates  the  estimation  of  the  fractal  dimension  for  most  of  the  liver 
ROIs. 


For  the  estimation  of  the  PSM,  the  range  of  the  values 
of  the  parameter  k ,  over  which  the  slope  of  the  best  fitting 
line  is  computed,  was  [15,20].  Additionally,  the  ranges  of 
the  values  of  the  parameter  k  used  for  the  estimation  of  the 
BCM,  MFE  and  K-NN  were  [7,19],  [3,15]  and  [50,150], 
respectively.  It  must  be  pointed  out  that  these  ranges  were 
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finally  selected  as  the  optimum  ones  in  relation  to  the  fractal 
dimensions  estimated  by  the  methods. 

From  Table  I,  it  can  be  observed  that  the  K-NN  may 
provide  an  indication  of  the  type  of  the  hepatic  lesion  under 
study,  based  on  its  fractal  dimension.  The  values  of  the 
fractal  dimension  using  the  other  two  methods,  BCM  and 
MFE,  are  quite  similar  and  they  could  not  provide  any 
information  for  the  discrimination  of  these  specific  hepatic 
lesion.  Nevertheless,  the  PSM  method  can  not  provide  any 
discrimination  between  the  two  categories  since  the  fractal 
dimensions  obtained  are  above  the  value  of  3;  thus  ROIs 
have  been  eliminated  from  the  study.  For  these  reasons,  a 
discriminant  analysis  is  then  performed. 

Furthermore,  a  well-known  algorithm  namely  the  Fuzzy 
C-Means  algorithm  [10]  was  then  applied  for  the  clustering 
of  the  input  data  into  two  clusters.  The  FCM  algorithm 
minimises  the  sum  of  squared  distances  to  the  prototypes 
weighted  by  constrained  membership  that  can  be  interpreted 
as  degrees  of  sharing.  The  function  has  two  arguments.  The 
first  is  a  feature  vector,  which  is  formed  by  all  the  values  of 
the  fractal  dimension,  to  be  classified  into  two  clusters.  The 
number  of  clusters,  which  is  the  second  argument, 
corresponds  to  the  normal  parenchyma  and  the 
hepatocellular  carcinoma,  respectively.  It  returns  the  cluster 
number  for  which  the  pattern  vector  has  the  highest  grade. 
The  FCM  clustering  for  each  hepatic  lesion  and  for  each 
fractal  dimension  estimator  is  then  presented  in  Table  II. 

Table n 

FCM  CLUSTERING  FOR  EACH  HEPATIC  LESION  AND  FOR  EACH 

_ FRACTAL  DIMENSION  ESTIMATOR _ 

BCM  MFE  K-NN 


Clusters 

Normal 

HCC 

Normal 

HCC 

Normal 

HCC 

Cluster  I 

48/99 

17/50 

72/99 

37/50 

92/99 

9/50 

Cluster  II 

51/99 

33/50 

27/99 

13/50 

7/99 

41/50 

Success 

rate 

48.48% 

66.00% 

72.73% 

26.00% 

92.93% 

82.00% 

From  the  results  of  the  Table  II,  it  is  obvious  that  the  k- 
th  nearest  neighbour  (K-NN)  outperforms  the  BCM  and  the 
MFE  methods.  Specifically,  the  K-NN  method  discriminates 
almost  93%  of  the  normal  parenchyma  (92  out  of  99  normal 
parenchyma)  and  82%  of  the  hepatocellular  carcinoma  (41 
out  of  50  hepatocellular  carcinoma)  correctly  to  the  two 
clusters.  Therefore,  the  K-NN  method  may  be  used  to 
discriminate  the  normal  parenchyma  from  the  hepatocellular 
carcinoma  using  CT  liver  images 

V.  Conclusion 

In  this  paper,  a  comparative  study  is  performed  to 
discriminate  two  different  CT  liver  lesions:  normal 
parenchyma  and  hepatocellular  carcinoma.  The  fractal 
dimension  of  99  ROIs  of  normal  images  and  of  50  ROIs  of 
hepatocellular  carcinoma  had  been  computed  throughout 
four  different  estimators,  regarded  as  the  most 
representative:  the  Power  Spectrum  Method  (PSM),  the 


Box-Counting  Method  (BCM),  the  Morphological  Fractal 
Estimator  (MFE),  and  the  kth-Nearest  Neighbour  estimator 
(K-NN),  as  a  novel  implemented  method.  The  results  of  the 
discriminate  analysis  show  that  the  K-NN  could 
discriminate  the  two  liver  lesions  more  efficiently  than  the 
BCM  and  MFE  methods.  In  particular,  the  PSM  was  proved 
to  be  insufficient  for  the  current  study,  since  underestimates 
the  fractal  dimension  of  most  of  the  CT  data.  The  success 
rate  of  discrimination  with  the  K-NN  method  reached  93% 
for  the  normal  parenchyma  and  82%  for  the  hepatocellular 
carcinoma.  In  the  future,  the  proposed  method  might  be  used 
for  the  extensive  analysis  of  other  liver  lesions  or  organs 
towards  to  the  classification  problem. 
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